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The general expression with the physical significance and positive definite 
condition of the eigenvalues of 4 x 4 Hermitian and trace-one matrix are 
, obtained. This implies that the eigenvalue problem of the 4x4 density matrix 

m . 
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, for an arbitrary state of two qubits is then given out and it is very easy 

o 

to use. Furthermore, we discuss some applications to the calculation of the 

Q-i! entanglement, the upper bound of the entanglement, and a model of the 
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transfer of entanglement in a qubit chain through a noisy channel. 
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is generally solved. The obvious expression of Peres' separability condition 
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The density matrix (DM) was introduced by J. von Neumann to describe the statistical 
concepts in quantum mechanics |Q. The main virtue of DM is its analytical power in 
the construction of the general formulas and in the proof of the general theorems. The 
evaluation of averages and probabilities of the physical quantities characterizing a given 
system is extremely cumbersome without the use of density matrix techniques. Recently, 
the application of DM has been gaining more and more importance in the many fields 
of physics. For example, in the quantum information and quantum computing H], DM 
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techniques have become an important tool for describing and characterizing the measure of 
entanglement, purification of entanglement and encoding [0-0]. However, even if DM is a 
simple enough 4x4 dimensional one, to write a general expression of its eigenvalues in a 
compact form with the physical significance seems not to be a trivial problem. Although one 
has known the theory of quartic equation, this is still difficult since DM has 15 independent 
parameters. Actually, we need a physical closed form for it but not a mathematical closed 
form only with formal meaning. This letter is just devoted to this fundamental problem in 
quantum mechanics. It successfully finds out the general expression of the eigenvalues of 
a 4 x 4 density matrix with a clear physical significance and in a compact form. Thus, an 
obvious expression of Peres' separability condition is derived clearly. This provides a very 
easy and direct way to use it. Moreover, some important applications to the entanglement 
and separability in the quantum information, such as the calculation of the entanglement, 
the upper bound of the entanglement and the transfer of the entanglement, are discussed 
constructively. 

The elementary unit of quantum information is so-called "qubit" [|7| . A single qubit can 
be envisaged as a two-state quantum system such as a spin-half or a two-level atom. A pair 
of qubits forms the simplest quantum register which can be expressed by a 4 x 4 density 
matrix. Just is well known, the eigenvalues of DM of two qubits are closely related with its 
entanglement and separability. For example, Wootters gave a measure of the entanglement in 
terms of the eigenvalues ||, and Peres' separability condition depends on the positive definite 
property of the partial transpose of DM. Therefore, it is very interesting and essentially 
important to know what is the general expression of eigenvalues of DM of two qubits in an 
arbitrary state. 

DM of two qubits can be written as 



where <Jq is two dimensional identity matrix and (7j is the usual Pauli matrix, p = p> 
(Hermitian) leads to a^ u be the real numbers, Trp = 1 (trace-one) requires aoo = 1, and 
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(1) 



from the eigenvalue of Pauli matrix it follows that — 1 < a^ v < 1. Moreover, it is easy to get 

= Tr(pcr M ® a v ). (2) 

Note that Eq.(|l]) does not involve with the positive definite condition for p. In order to 
find out the general expression of the eigenvalues of DM, we first give out the following two 
lemmas. 

Lemma One The form of the characteristic polynomial of a 4 x 4 Hermitian and trace-one 
matrix Q is 

b + M + M 2 - A 3 + A 4 , (3) 
where the coefficients b , bi and 6 2 are defined by 

bo = ^[1 ~ t^B ~ (M A ) 2 ~ (Mb? 
+2&A£b+ ((^ A ) 2 -TrA 2 )£ A - £ B 
+2f B A 2 $, A - 2TrA f B A£ A - (a x x a 2 ) 2 
-(o a x a 3 ) 2 - (o 3 x a x ) 2 - 2(oi x a 2 ) • a 3 ] 

-1[™ 2 -(M 2 ) 2 ], (4) 
&! = i[2Trfi 2 - 1 - f A A£ B + ( ai x a 2 ) • a 3 ], (5) 

o 

&2 = ^(1-M 2 ). (6) 

In the above equations, we have introduced the polarized vectors of the reduced density 
matrices $, A = (a w , a 2 o, 030), £b = ( a oi, a 02, 003); the space Bloch's vector a,\ = (an, a i2 , 013), 
a 2 = (a 2 i, a 22 , a 23 ), a 3 = (031,032,033); and the polarized rotation matrix A = {a^} = 
1, 2, 3). Note that £{a,b} * s viewed as a column vector and its transpose $Ja,b} * s then a row 
vector. The physical meaning of 3 x 3 matrix A can be seen in my paper H. Again, the 
positive definite condition has not been used here and a^ v is defined just as Eq.(|2]). 

Lemma Two If a 4 x 4 Hermitian and trace-one matrix Q has m non-zero eigenvalues, then 

TrQ 2 >-. (7) 

m 



If Q is positive definite, then 



TrQ 2 < 1. 



(8) 



To prove Lemma one, we need to use the physical ideas to arrange those coefficients from 
the characteristic determinant into a compact form. And, it leads to the general expression 
of the eigenvalues of DM with the physical significance for their applications. Lemma two 
can be obtained by the standard method to find the extremum. Based on the theory of the 
quartic equation, we have 

Theorem One The eigenvalues of the 4x4 Hermitian and trace-one matrix Q are 
1 1 



^(-) = 7 - 



4 AV3 
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2^6 
1 1 



(4Trft 2 - 1 + 8ci cos 0) 1/2 

3^3(1 + 8&i - 2TrQ 2 ) 



4TrST -1 -Act cos0 + 



4 4V3 
1 



v/4Trfi 2 - 1 + 8ci cos ( 
(4M 2 - 1 + 8ci cos 0) 1/2 



1/2 



(9) 



4Trfi 2 - 1 - 4ci 



COS ( 



3^(1 + 86i - 2TrQ 2 ) 

OCi COS( 



1 1/2 



(10) 



where 



cos . 



C 2 a 3 

25 = 4C08 ' 



3 cos ( 



126 + 36i + b 2 



2- 



c 2 = 27b\ + 6 (27 - 72& 2 ) + 9&i& 2 + 2& 3 . 



(11) 
(12) 
(13) 



In the above equations, we have assumed c\ ^ 0, c 2 ^ 0. Note that c 2 — 4c? = — 27(Ai _ 
A 2 ) 2 (Ai - A 3 ) 2 (A! - A 4 ) 2 (A 2 - A 3 ) 2 (A 2 - A 4 ) 2 (A 3 - A 4 ) 2 is non-positive. Thus, 4c? > c 2 > 0, 
and so c\ is real. If there is any repeated root (c 2 , — 4c? = 0), <fi = or 7r/3 since c 2 = ±2cf . 
In fact, Eq.(p|) implies that 



cos< 



2 2 /3(c 2 + v/c 2 - 4cf)V3 



+ 



(02 + y^^f) 1 / 3 

2 x 2V3 Cl 



Arg 



c 2 + 



4c? 



1/3" 



(14) 
(15) 



Obviously, if c 2 < 0, then n/6 < < ir/3, and if c 2 > 0, then < < n/6. 

Now consider the case that C\ or/and c 2 are equal to zero. We discuss it by two steps. 
First, suppose b 2 = 3/8 or Trf2 2 = 1/4. If only c\ = or c 2 = 0, then some of these 
eigenvalues will be complex numbers. This is contradict with Hermitian property of DM. 
So this conclusion means that both c\ and c 2 have to be zero together. Thus, we can obtain 
that all the eigenvalues are 1/4. Second, suppose 6 2 ^ 3/8 or Trfi 2 ^ 1/4. We have to 
analysis the possibilities stated as the following. 

For only a = or 12b + 3&i + b\ = 0, again from c| = c 2 . - 4c? = -27(Ai - A 2 ) 2 (Ai - 
As) 2 (Ai — Azt) 2 (A 2 — As) 2 (A 2 — A4) 2 (A3 — A4) 2 < 0, we have that c 2 has to be zero since it is 
real. So only c\ = is impossible. 

For only c 2 = or 27bj + b (27 - 72b 2 ) + 96i& 2 + 26| = 0, the eigenvalues of the 4 x 4 
Hermitian and trace-one matrix are then 

/ r- \ !/ 2 

'*( ) 1 \ \(— 1 1 3\/3(l + 8^1 — 2Tr^ )\ 
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l v '3 ^ "" V4Trfi 2 - 1 

A±(+ ) = I + J-v^TIiJ- C 4T,tf - 1 - 3^(1 +8^-2™ 
4 4^ 2V6 \ V4Trfi 2 - 1 



2, x V 2 



(17) 



For both c\ — and c 2 = 0, the eigenvalues of the 4x4 Hermitian and trace-one matrix 
are either 

Aias = \- ^V4W - 1, (18) 

A 4 = - + ^V4M 2 - 1. (19) 
4 4 

if 60 = [3-6Trfi 2 -6(™ 2 ) 2 + v / 3(4M 2 -l) 3 / 2 ]/288, 61 = [18M 2 -9- v / 3(4™ 2 -l) 3 / 2 ]/72, 
or 

Ai, 2 ,3 = \ + ^=V4W - 1, (20) 

A 4 = - - ^V4M 2 - 1. (21) 
4 4 

if 60 = [3-6™ 2 -6(™ 2 ) 2 - v / 3(4™ 2 -l) 3 / 2 ]/288, 61 = [18Trfi 2 -9+V3(4Trft 2 -l) 3 / 2 ]/72. 
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Just is well known, Peres' separability condition tells us, all the eigenvalues of the partial 
transpose of DM ought to be non- negative ||. Thus, taking the minimum eigenvalue in 
Theorem one and setting it non-negative, we have 

Theorem Two The separability condition of DM p of two qubits in an arbitrary state is 

1 > -^(4Trp 2 - 1 + 8cf cos <f) 1/2 + ^= [4Trp 2 - 1 
V3 v 6 



. P , P 3^(1 + 86? - 2Trp 2 ) 
— 4q cose 



V^Trp 2 - 1 + 8cf cos0 p 



(22) 



where 



&o = bo ~ ^[((TrA) 2 - TrA 2 )£ A • £ B + 2^A 2 ^ 

-2TrA + i( 0l x oa) ■ a 3 , (23) 

6f = 6i - -(oi x a 2 ) ■ a 3 , (24) 

62 = 62, (25) 



c? = ^126^ + 36^ + 6^ (26) 



4 = 276f 2 + 6^(27 - 726^) + %{ b v 2 + 26™ (27) 
cos0 = = ,„ p . (28) 

22/3(^ + ^2_ 4c P6)l/3 2X2V3 C P ^ > 



And q ^0, ^0. If only c£ = 0, the separability condition becomes 



1/2 



If both cf = and cf = 0, then in case one the DM is always separable and in case two the 
separability condition is 

Trp 2 < i (30) 

In the above, we have used the fact that the trace of the square of the partial transpose 
matrix of DM is equal to the trace of the square of DM. 

Obviously, the pure state is the simplest case. In fact, we can prove the following theorem: 
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Theorem Three The eigenvalues of the partial transpose of DM of two qubits in a pure state 
|0) = a |oo) + 6|01) + c|10) + is 

=F \ad - bc\, hi =F y/l - A\ad - bc\ 2 ), (31) 

and then the separability condition is just 

ad-bc = 0. (32) 

It is consistent with my paper ||. 

Because Peres' separability condition is necessary and sufficient one for two qubits, the 
Theorem two and Theorem three, as the obvious and general expression of Peres' condition, 
are necessary and sufficient one either. 

If there are some vanishing eigenvalues for a 4 x 4 Hermitian and trace-one matrix, the 
conclusions can be simplified. The following theorems will show this judgment. Their proofs 
can be given by solving the corresponding characteristic equations. 

Theorem Four If at least there is one vanishing eigenvalue for a 4 x 4 Hermitian and trace-one 
matrix, its eigenvalues are 

Ai = hi + a/6M 2 - 2 cos0), (33) 
A 2 = -[1 - V6TrQ 2 - 2 cos(0 - tt/3)], (34) 
A 3 = -[1 - VGTrQ 2 - 2 cos(0 + tt/3)], (35) 



where 



VT^3&2 

COSffl 



2 2 / 3 (rf+ v /d 2 -4(l-36 2 ) 3 )V3 
(d + v /rf 2 -4(l-36 2 ) 3 ) 1 / 3 



(36) 



2 x 2 1 /3 v /i - 36 2 
d = 2- 276i - 96 2 . (37) 

Here we have assumed 3TrQ 2 — 1 ^ and d = (3Ai — 1)(3A 2 — 1)(3A3 — 1)^0. If d < 0, then 
vr/6 < (f) < tt/3, and if d > 0, then < (fx tt/6. Because that d 2 - 4(1 - 36 2 ) 3 = -27(A X - 



A 2 ) 2 (A 2 - A 3 ) 2 (A 3 - AO 2 , we have 4(1 - 36 2 ) 3 > d 2 > 0, and 1 - 36 2 = (3M 2 - l)/2 > 0. If 
Trf2 2 = 1/3, we have that d has to be zero. Thus, &i = —1/27 and b 2 = 1/3. This implies 
that all the eigenvalues are equal to 1/3. In particular, only d — 0, the eigenvalues becomes 

Ai = \, (38) 
A 2 , 3 = ^1 ± y|v3 W - . (39) 

Theorem Five If at least there is one vanishing eigenvalue for a 4 x 4 Hermitian and trace-one 
matrix, then the positive definite condition of eigenvalues is 

V6M 2 - 2 cos(0 p - tt/3) < 1. (40) 

If only d = 0, the positive definite condition becomes 

M 2 < jj. (41) 

Theorem SixU at least there are two vanishing eigenvalues for a 4x4 Hermitian and trace-one 
matrix, then the other eigenvalues are 

1 



A± = -(1 ± V2M 2 - 1). (42) 
2 

The positive definite condition of the eigenvalues 

Trf2 2 < 1. (43) 

Thus, from the Peres' separability condition, it is easy to prove the following theorem: 

Theorem Seven If the partial transpose of DM of two qubits has at least two vanishing 
eigenvalues, this density matrix is separable. If the partial transpose of DM of two qubits 
has only one vanishing eigenvalue, the separability condition is obtained by using of Theorem 
Five to it and setting the minimum eigenvalue is positive. 

Now, let's we discuss some applications of our theorems. Just is well known that many 
measures of entanglement are related with the quantum entropies which are defined by 
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density matrix, for example, the entanglement of formation |10| and the relative entropy 
of entanglement JO]] . To compute quantum entropy, we often need to find the eigenvalues 
of density matrix. Even, according to Wootters ||, the measure of entanglement of two 
qubits is directly determined by the eigenvalues of 4 x 4 hermitian matrix. Therefore, in 
terms of our Theorems about the eigenvalues, we can easily calculate the entanglement of 
formation of an arbitrary state of two qubits. As to the relative entropy of entanglement or 
its improving ||12j| , we have to calculate von Neumann entropy which is just defined directly 
by the eigenvalues. 

Furthermore, we can find a relation between the upper bound of entanglement and the 
eigenvalues. 

Theorem Eight If all the eigenvalues of DM of two qubits are not zero, the possibly maximum 
value of the entanglement of formation is not larger than 

-^(4Trp 2 - 1 + 8cf cos P ) 1/2 + -= [4Trp 2 - 1 
V3 v 6 



P P 3^(1 + 86f - 2Trp 2 ) 



-4c[ cos + 



1 1/2 

(44) 



v/4Trp 2 - 1 + 8cf cos0 p 
This is because DM of two qubits can be written as 

P = KiJ + (1 - 4A min )p'. (45) 

Note that we can not put a number larger than A m i n in front of the identity matrix J, because 
we have to keep p' to be positive definite. 

Furthermore, let's consider a model of transfer of entanglement jDJ. This can be ex- 
pressed as a following story. Alice and Bob are friends. One day, Alice and Bob sat together 
at the lounge in a party. On the left side of Alice is Charlie and on the right side of Bob 
is David. Alice and Charlie, Bob and David respectively exchanged their seats. This leads 
to Alice and Bob's entanglement decreases. In language of quantum information, Alice and 
Bob shares an entangled state initially. Without loss of generality, suppose it is in Bell state 
(|00) + \11))/V2, Charlie is in |c) and David is in \d). That is that four of them is in a total 
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state |c) g) ( 1 00) + |11)) £g> \d)/y2. Now, introduce the swapping interaction [|T4| respectively 
between Alice and Charlie, and between Bob and David: 



/ 



S 



10 
10 
10 
1 



\ 



S\ab) = \ba) (a, b = 0,1). 



(46) 



It is easy to see that 
1 



S® S 



V2 



c)<g>(|00) + |ll))® \d) 



V2 



(|0) ® \cd) |0) + |1) g> |cd) <g> |1)). (47) 



Thus, the entanglement between the second qubit and the third qubit is transfered to the 
entanglement of the first qubit and the fourth qubit. In general, the swapping process suffers 
the affection of noisy. We suppose, after transfer of entanglement, that DM becomes 

1 



P 



e)p 



e 4 J - 



(48) 



where e represents the strength of the noise and p is, in form, the same as DM before transfer 
of entanglement. Obviously, this model can be extended to a qubit chain: 



Pn 



pi 



"V" 



n 



(49) 



Through the swapping interaction, the entanglement can be transfered along with the chain 
one node by one node and forward to the opposite directions. At the beginning, denote DM 
for a pair of given adjacent nodes of a qubit chain as po- After the first swapping, the DM 
of a pair of qubits respectively on the nodes of the left side and right side of the given pair 
of qubits is written as p\. Since the affection of noise, after the n-th swapping in turn along 
with two directions, the DM of a pair of qubits respectively on the n-th nodes of the left 
side and right side of the original two adjacent qubits becomes 



p n = (1 - e)p n -i + e-I. 



(50) 
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This equation has not taken into account the fact that p n and p n _i are related with the 
different pairs of qubits, and let it is valid only in mathematics. We also assume that at 
the beginning a pair of qubits in the adjacent nodes is in a pure state. Thus we would like 
to know what is n's value if p n is separable. That is to calculate the transfer distance n of 
entanglement for such a chain of qubits through a noisy channel. 

According to our theorem, the minimum eigenvalue of the partial transpose of p n is 

A min = - ((1 - e) n (l + 4\ad - bc\)}. (51) 

By using of Peres' criterion for the separable state, we can find that 

n < _'°8(1 + 4M-M) (52) 
log(l - e) 

In particular, when p is a density matrix in the maximum entangled state, we obtain 

n < — - ; r. 53 

- log(l-e) 1 J 

Obviously, when one hopes n = 10, then it allows the noisy strength e is not larger than 
0.104042. When the noisy strength e is larger than 0.42265, any transfer will lead in disen- 
tanglement. If the noisy strength e only reaches at 0.01 or 0.1, the transfer distance n can 
be 109 or 10. Likewise, we can describe the transfer of entanglement along with one direc- 
tion. The significances and applications of this model of transfer of entanglement should be 
imaginable. 

In addition, we hope to apply our theorems into seeking the minimum pure state decom- 
position of DM, this work is in progressing. In a words, we can say, of course, the theorems 
proposed here are the useful tools to study the entanglement and the related problems. 

I would like to thank Artur Ekert for his great help and for his hosting my visit to center 
of quantum computing in Oxford University. 
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